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1.  MODELING  THE  REFRACTIVE  INDEX  STRUCTURE  PARAMETER 
IN  THE  MARINE  PLANETARY  BOUNDARY  LAYER 


1.1  INTRODUCTION 

At  the  request  and  under  the  sponsorship  of  the  Navy's  High 

Energy  Laser  (HEL)  Program  Office,  the  Naval  Environmental 

Prediction  Research  Facility  (NEPRF)  hosted  a Workshop  on  Modeling 

2 

of  Refractive  Index  Structure  Parameter,  in  the  Marine 

Planetary  Boundary  Layer,  Appendices  A and  B are  list  of  attendees 
and  the  agenda  at  the  workshop. 

2 

NEPRF  proposed  a preliminary  model  to  serve  as  a strawman 

for  discussion  at  the  workshop.  The  model  is  included  as  Appendix 
C.  The  workshop  participants  reviewed  the  components  of  the 
preliminary  model,  and  suggested  modifications  to  improve  the 
mo  del  . 

2 

This  document  summarizes  the  modified  C model,  and  presents 

n 2 

a state  of  the  art  review  of  marine  planetary  boundary  layer 
modeling.  The  document  represents  the  consensus  opinion  of  the 
workshop  participants. 

1 .2  HEL  PROBLEM 

1.2.1  Atmospheric  Effects  on  HEL  System  Performance 

Atmospheric  conditions  degrade  HEL  system  performance  by 
reducing  fl uence  on  target  in  a variety  of  ways;  aerosols  and 
water  vapor  absorb  energy;  aerosols  scatter  energy  out  of  the  beam 
p.’th;  atmospheric  turbulence  spreads  the  beam  fluence;  and  atmos- 
pheric turbulence  also  causes  the  beam  to  wander  off  its  intended 
target.  Figure  1.1  is  a schematic  representation  of  these  four 
effects  which  are  shown  separately  and  sequentially  in  the  figure. 
In  nature,  each  effect  occurs  simultaneously  along  the  beam. 

The  aerosol  and  water  vapor  elements  of  the  problem  are  being 
considered  elsewhere  within  the  Navy  R&D  commi^nity.  Consequently, 
these  elements  were  not  addressed  at  the  NEPRF  workshop,  nor  do 
we  consider  them  in  this  review. 


occur  together  and  anywhere  along  the  beam. 


1.2.2  Turbulence  Near  the  Sea  Surface 

The  Navy  is  developing  laser  systems  for  shipboard  operation. 
From  a ship,  a laser  could  be  fired  horizontally,  vertically  or  at 
any  angle.  In  the  first  case,  turbulence  near  the  sea  surface 
degrades  the  beam.  In  the  latter  two  cases,  turbulence  aloft  as 
well  as  near  the  surface  can  degrade  performance. 

The  degradation  of  system  performance  is  caused  by  turbulent 
fluctuations  of  refractive  index,  weighed  nonlinearly  toward  the 
transmitting  end  of  the  propagation  path.  Therefore  turbulent 
conditions  affect  beam  quality  more  when  they  are  near  the  source 
than  similar  conditions  near  the  target.  We  call  this  range 
dependence.  System  performance  calculations  take  range  dependence 
into  consideration. 

For  a surface  based  laser  system,  range  dependence  makes 
turbulence  near  the  ground  degrade  system  performance  more  than 
turbulence  aloft.  For  this  reason  and  others  to  be  elaborated 
upon  later,  we  focus  upon  the  atmosphere  near  the  sea  surface.  We 
call  this  region  the  marine  planetary  boundary  layer  (MPBL) 

(Fig.  1.2  and  1.3). 

1 .3  TECHNICAL  PROBLEM 

There  are  two  elements  of  the  technical  problem  that  must  be 
recognized  by  the  reader  of  this  report;  (1)  the  turbulence  param- 
eter that  is  related  to  HEL  system  performance  is  the  refractive 

2 2 

index  structure  parameter,  C^;  and  (2)  C^,  like  many  meteorological 
parameters,  behaves  differently  near  the  surface,  in  the  MPBL, 
than  it  does  aloft  in  the  free  troposphere. 

2 

1.3.1  Refractive  Index  Structure  Parameter,  C^ 

The  structure  parameter  for  any  atmospheric  variable  is  a 

measure  of  the  turbulent  fluctuations  of  that  variable  over  small 

distances.  For  HEL  systems,  the  atmospheric  variable  of  interest 

2 2 

is  refractive  index,  and  the  structure  parameter  is  C„.  C is 

n n 

defined  by: 

C^  = [n(x  + r)  - n(x)]^  r’^^^  (1  .1) 


3 


TEMPERATURE 


FREE  ATMOSPHERE 


MARINE  INVERSION 


Fig.  1.2  Schematic  Marine  Atmosphere  Temperature  and  Moisture 
Structure.  The  vertical  scale  is  non-linear  to  allow  more 
detail  in  the  marine  planetary  boundary  layer  (MPBL).  The 
MPBL  is  typically  approximately  1 km  thick,  well  mixed  as 
evidenced  by  constant  mixing  ratio  and  capped  by  a stable 
marine  inversion  at  a height  Z. . The  free  atmosphere  extends 
from  the  marine  inversion  to  the  tropopause  and  is  charac- 
terized by  minimal  turbulent  mixing. 


FREE  ATMOSPHERE 


UPPER  MPBL 
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SURFACE  LAYER 


Fig.  1.3  Marine  Planetary  Boundary  Layer  (MPBL).  The  surface 
layer  is  approximately  10%  of  the  MPBL  and  is  a region  of 
surface  forced  mixing.  The  upper  MPBL  is  influenced  by  both 
surface  forcing  and  free  atmospheric  effects.  Z.  is  the 
height  of  the  marine  inversion  that  caps  the  MPBil  . 


5 


n(x+r)  and  n(x)  in  (1.1)  are  refractive  index  at  points  in  space 

sepa rated  by  a short  distance,  r,  and  the  overbar  represents  an 

average  over  enough  time  and/or  space  to  obtain  a statistically 

stable  sample  of  fluctuation.  In  the  atmosphere,  the  averaging 

time  for  such  a sample  is  approximately  30  minutes  for  a fixed 

point  measurement. 

2 

then  describes  the  intensity  of  refractive  index  fluctua- 
tions over  a range  of  the  size  of  turbulent  eddies.  The  eddies 
affecting  laser  propagation  most  are  determined  by  laser  wave- 
lengths, beam  diameter,  and  beam  range.  For  the  present  candidate 
system,  the  eddy  size  is  roughly  10  to  20  cm. 

Refractive  index,  n,  is  a secondary  atmospheric  variable 
that  is  derived  from  a combination  of  the  primary  variables 
pressure,  P,  temperature,  T,  and  water  vapor  pressure,  e.  The 
exact  combination  of  the  primary  variables  required  to  compute 
refractive  index  depends  upon  the  wavelength  of  the  energy  being 
refracted.  In  the  case  of  a laser  or  any  other  optical  system, 
the  equation  for  refractive  index  is  by 

n-1  = + A2f  (1  .2) 

where 

A^  = 78.7  X lO'^K  mb"^  and  A^  = 66.3  X lO'^K  mb'^  . 

Since  refractive  index  is  related  to  temperature  and  moisture 

2 

the  fluctuations  in  refractive  index  that  compose  C are  similarly 
composed  of  fluctuations  in  temperature  moisture  and  pressure. 

The  effect  of  pressure  fluctuations  can  be  neglected.  The 

fluctuations  of  temperature  and  moisture  are  described  by  tempera- 

2 . 2 
ture  structure  parameter,  C^,  moisture  structure  parameter,  C^, 

and  a structure  parameter  relating  the  two  called  temperature/ 

moisture  structure  parameter,  '^eT- 

The  relationship  between  the  individual  component  structure 
2 

parameter  and  C is 


where  , B^.  and  B^  are  functions  of  temperature  pressure,  vapor 
pressure . 

Over  the  range  of  normally  observed  values  of  the  various 

2 ? 
structure  parameters,  Cy  is  the  largest  contributor  to  C^.  Over 

land,  where  moisture  sources  are  more  remote  and  relatively  small, 

2 

the  contributions  by  and  C^j  are  negligible. 

Most  observational  studies  of  the  structure  parameter  have 

been  made  over  land.  Consequently,  of  the  three  structure  param- 

2 

eters,  we  know  the  most  about  C^-.  Because  the  contribution  of  the 

2 

other  structure  parameters  is  negligible  over  land  anyway,  C is 

2 ^ 
solely  determined  by  Cy  over  land  (e.g.,  over  White  Sands  Missile 

Range)  . 

Over  water,  however,  the  contributions  of  the  moisture 

2 

structure  parameters  are  significant.  When  we  desire  C , we  can 

2 

observe  and/or  calculate  Cy  and  then  adjust  it  for  the  contribu- 
tion made  by  fluctuations  in  moisture. 

1.3.2  Atmospheric  Layers 

The  marine  atmosphere  is  composed  of  two  layers;  the  marine 
planetary  boundary  layer,  previously  defined  as  the  MPBL,  and  the 
free  atmosphere.  The  MPBL  extends  from  the  ocean  surface  to  a 
stable  region  typically  about  1 km  above  the  surface  that  we  call 
the  marine  inversion.  The  free  atmosphere  extends  from  the  marine 
inversion  to  the  tropopause,  another  stable  region  at  a height  of 
approximately  15  km.  Figure  1.2  shows  a schematic  marine  atmos- 
pheric temperature  and  moisture  sounding  showing  the  MPBL,  marine 
inversion,  free  troposphere  and  tropopause. 

1.3. 2.1  Free  troposphere 

Motion  in  the  free  atmosphere  is  controlled  by  large  scale 
temperature  gradients.  On  the  global  scale,  air  rises  near  the 
equator,  cools  as  it  moves  north  and  descends  in  the  mid  latitudes. 
On  a synoptic  scale  in  the  Northern  Hemisphere,  air  moves  clockwise 
around  and  descends  in  high  pressure  systems,  and  moves  counter- 
clockwise around  and  ascends  in  low  pressure  centers. 


Surface  effects  exert  little  or  not  influence  on  motion  in  the 
free  troposphere.  With  the  exception  of  clear  air  turbulence 
(CAT),  there  is  little  atmospheric  turbulence  in  the  free  tropo- 
sphere. Turbulent  fluctuations  of  re fr active  index,  are  weak  in 

2 

the  free  troposphere  and  consequently  values  are  low  in  this 
region,  except  possibly  near  CAT.  Since  CAT  accompanies  strong 
wind  shears  and  sharp  temperature  discontinuities,  CAT  occurs  less 
than  20%  of  the  time. 

Potential  Navy  HEL  weapons  are  intended  to  operate  from 
ships.  Since  turbulence  effects  on  HEL  systems  performance  are 
range  dependent,  MPBL  turbulence  is  the  primary  turbulence  influence 
on  laser  propagation,  while  effects  of  the  free  atmosphere  are 
minimal.  For  this  reason,  the  NEPRF  workshop  considered  only  MPBL 
turbulence.  The  workshop  did  not  consider  free  atmosphere  turbu- 
lence nor  will  we  review  it  in  this  report. 

We  suggest,  however,  that  a search  be  made  of  the  literature 
to  confirm  the  conditions  for  low  values  of  Cj^  in  the  free  tropo- 
sphere and  that  HEL  system  performance  calculations  of  the  effects 

2 

of  different  levels  of  C^  at  free  troposphere  distances  from  the 
weapon  be  made  to  confirm  the  basic  premises  upon  which  we  excluded 
the  free  troposphere  from  this  report. 

1.3. 2. 2 Marine  planetary  boundary  layer 

In  contrast  to  the  free  troposphere,  motion  in  the  MPBL  is 
only  weakly  influenced  by  large  scale  temperature  discontinuities 
and  is  heavily  influenced  by  surface  effects.  The  primary  surface 
effects  are  heating  or  cooling  and  friction. 

Surface  effects  cause  the  PBL  to  be  highly  turbulent. 

Surface  heating  causes  buoyancy  generated  turbulence  and  wind 
shear  causes  mechanically  generated  turbulence. 

The  MPBL  can  be  divided  into  a layer  in  immediate  contact 
with  the  surface,  called  the  surface  layer,  and  the  layer  above 
that  up  to  the  marine  inversion,  called  the  upper  MPBL.  The 
surface  layer  comprises  approximately  \0%  of  the  ^PBL  and  experi- 
ences almost  no  free  tropospheric  influence.  The  upper  MPBL  can 
be  thought  of  as  a transition  layer  between  the  surface  layer  and 
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the  free  atmosphere,  where  both  surface  forcing  and  synoptic  scale 

2 

features  influence  C^.  Figure  1.3  is  a schematic  representation 
of  the  marine  PBL  showing  the  relative  depth  of  the  surface  layer 
and  the  upper  MPBL  . 

1.3. 2. 2.1  Surface  layer.  Scaling  laws  relating  to  bulk 
meteorological  parameters  in  the  surface  layer  have  been  success- 
fully studied  for  a number  of  years  and  are  well  defined. 

Section  2 of  this  report  describes  surface  layer  scaling  laws 
and  describes  a sequence  of  calculations  to  compute  from  ship- 
board measurements  of  air  temperature,  relative  humidity,  wind 
speed  and  sea  surface  or  bucket  temperature.  It  was  the  consensus 
of  the  workshop  that  the  surface  layer  scaling  laws  predict  to  a 
factor  of  two.  Further  refinement  is  unjustified  because  the 

accuracy  of  the  shipboard  observations  introduces  the  same  order 

2 

of  uncertainty  into  C^. 

1.3. 2. 2. 2 Upper  MPBL . Scaling  laws  in  the  upper  MPBL  are  only 
qualitatively  known.  Many  of  the  relevant  physical  processes  that 
control  turbulent  fluctuations  in  the  upper  MPBL  are  understood; 
however,  several  are  not.  The  largest  unknown  influences  are 
those  that  feed  down  from  the  free  atmosphere. 

In  the  upper  MPBL,  two  avenues  of  approach  are  available  to 

link  C to  bulk  meteorological  variables:  (1)  scaling  laws  can 

be  developed  to  the  point  where  they  are  as  well  accepted  as  those 

in  the  surface  layer;  or  (2)  scaling  laws  can  be  abandoned  in 

2 

favor  of  dynamic  models  of  the  upper  MPBL  that  calculate 
directly  from  first  principles. 

Chapter  3 of  this  report  describes  the  current  status  of 
both  scaling  laws  and  dynamic  models  in  the  upper  MPBL.  It 
concludes  that  considerable  progress  needs  to  be  made  to  improve 
scaling  laws,  and  to  verify,  and  probably  improve,  dynamic  models 
in  the  upper  MPBL.  It  further  concludes  that  dynamic  models  may 
be  used  to  develop  improved  scaling  laws  in  lieu  of  large  numbers 
of  observations  . 
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2.  ESTIMATION  OF  REFRACTIVE  INDEX  STRUCTURE  PARAMETER  IN  THE 
SURFACE  LAYER  FROM  BULK  METEOROLOGICAL  MEASUREMENTS 

2.1  INTRODUCTION 

Modern  turbulence  theory  provides  a good  description  of 

fluctuations  of  wind,  temperature  and  humidity  in  the  atmospheric 

surface  layer,  i.e.,  the  layer  of  the  atmosphere  within  50  m to 

100  m of  the  surface.  With  this  theory,  the  relevant  turbulence 
2 

parameter,  C^,  is  related  to  commonly  observed  meteorological 
measurements,  such  as  temperature  and  wind,  by  well  accepted 
scaling  laws.  A preliminary  NEPRF  model  that  incorporated  these 
scaling  laws  in  a continental  environment  was  supplied  to  the 
workshop  and  is  included  as  Appendix  C.  The  following  paragraphs 
present  a summary  of  the  scaling  laws  and  how  they  are  adapted  to 
the  marine  environment. 

2 

Paragraphs  2.2  through  2.8  describe  the  relationships  of 

to  a progression  of  physical  scaling  laws  which  eventually  depend 

on  observable  meteorological  measurements.  Paragraph  2.9  describes 

the  modification  of  one  measured  parameter,  sea  surface  temperature 

from  the  bucket  measurement,  to  actual  sea  surface  temperatures. 

Paragraph  2.10  summarizes  the  values  of  empirical  constants  in  the 

calculation.  A method  of  calculation,  with  a flow  diagram,  are 

2 

included  in  Para.  2.11.  The  problem  of  i ntermi  ttency  of  is 
addressed  in  Para.  2.11.  Recommendations  for  the  modification  of 
the  NEPRF  preliminary  model  are  presented  in  Para.  2.12. 

2.2  REFRACTIVE  INDEX  STRUCTURE  PARAMETER,  C^ 

2 

Earlier  we  have  stated  that  Cn  depends  principally  on  the 

2 

temperature  structure  parameter,  Cj,  with  a modification  for 
humidity  fluctuations  and  the  correlation  of  humidity  and  tempera- 
ture fluctuations.  At  the  NEPRF  workshop,  two  alternative 

formulations  of  the  humidity  modification  were  presented.  The 

2 

first,  introduced  by  Wesely,  modifies  Cj  by  a function  involving 
moisture  and  the  correlation  of  moisture  and  temperature 


11 


fluctuations.  The  second  approach.  Introduced  by  Panofsky,  uses 


conventional  surface  layer  scaling  laws  to  derive  the  moisture 

2 

modification  of  C„. 

n 


Both  methods  lead  to  a modification  of  temperature  structure 

function  to  account  for  humidity.  When  temperature  and  moisture 

fluctuations  are  perfectly  correlated,  the  modifications  are 

identical.  If  the  correlation  is  not  perfect,  Wesely  indicates  a 

smaller  humidity  correction  than  Panofsky. 

The  importance  of  the  existence  of  two  derivations  is  in 

2 

their  independent  verification  of  the  maximum  change  of  Cj  to 
account  for  moisture.  Together  with  experimental  verification  by 
Friehe  ( 1 977  ) and  Davidson  et  al  . ( 1 978),  these  derivations  justify 


confidence  in  the  calculation  of  the  moisture  modification  to  C 


2 


Wesely  (1976)  has  proposed  an  approach  to  evaluating  the 

2 

influence  of  humidity  on  C . Using  the  relationship  of  the 


refractive  index  of  air  to  pressure,  temperature, and  vapor  pressure 

introduced  in  1.5.1,  he  expands  the  definition  of  refractive  index 

structure  parameter  in  terms  of  temperature  and  specific  humidity 

2 2 

structure  functions,  Cj  and  Cg , and  the  correlation  of  temperature 
and  vapor  pressure,  C, 


eT- 


These  two  expressions  are  combined  to  provide  C as  a func- 


tion of  the  requisite  structure  functions  and  covariances. 


= c; 


2.  2 _ ^ ^ 2r 

■(1-A2/At  )CgT‘ 

■(1-A2/A,  )CJ' 

\T“;  %’  “v  - ' ^ 2rg^ 

L P^T  J 

L PC,  J 

(2 


where  T is  mean  layer  temperature,  P is  mean  layer  pressure 

C 

= 78.7-10‘^K  mb‘\  A2  = 66.7-10"®K  mb"^  and  r^.^  = 

This  relation  becomes  simple  when  '“eT  " ^ ’ that  is  e,  T are 
perfectly  correlated. 

Until  r^j  is  better  related  to  atmospheric  parameters,  we 
recommend  using  r ^ equal  to  1 and  hence  Eq . (2.2). 


Panofsky  proposed  an  alternative  derivation.  Based  on 
surface  layer  scaling  laws,  the  derivation  does  not  explicitly  use 
the  temperature-humidity  correlation,  C^y,  but  ends  with  the  same 
expression  as  Eq.  (2.2).  The  derivation  is  shown  in  Appendix  D. 

2.3  TEMPERATURE  STRUCTURE  PARAMETER  Cy 

2 

The  temperature  structure  parameter  Ct  is  defined  in  a similar 

2 

way  to  the  definition  of  provided  in  1.5.1, 

Cy  = [T(x+r)  - T(x)]^  r’2/^  (2.3) 

where  T(x)  is  the  temperature  at  point  x,  and  the  remaining  nota- 
tion is  the  same  as  in  1.5.1.  In  the  surface  layer,  dimensional 

arguments  first  introduced  by  Monin  and  Obukhov  can  be  used  to 
2 

relate  Cy  to  a temperature  scale,  T,,,  a height  z,  and  a stability 
function  g(z/L) , 

Cy  = T*^  z-2/3  g(z/L)  (2.4) 

where  L is  the  Moni n-Obukhov  length  scale. 


2.4  STABILITY  FUNCTION  IN  Cy 

Using  a set  of  careful  observations  over  uniform  terrain  in 
Kansas,  Wyngaard  et  al  . (1971)  have  empirically  determined  the 
function  g(z/L)  from  measurements  of  Cy  and  the  scaling  parameters 
L and  T* . The  results  have  been  tested  over  water  by  Friehe  (1977) 
and  Davidson  et  al  (1978).  As  a result,  the  following  expressions 
are  recommended: 


g(z/L)  = 4 .9(1-7z/L)’^^^  (L<0) 

g(z/L)  = 4.9(l+2.75(z/L)‘')  (L>0) 

where  the  exponent  b has  been  suggested  either  as  1 or  .667. 


(2.5) 
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2.5  SCALING  FUNCTIONS  IN  THE  SURFACE  LAYER 


The  Mont n-Obukhov  theory  describes  the  surface  layer  using 
non-dimensional  atmospheric  variables.  Physical  variables  of 
length,  z,  wind  speed  U,  temperature  T,  and  specific  humidity  q 
are  made  non-dimensional  by  defining  scaling  length,  L,  scaling 
speed  U*,  a scaling  temperature  T*  and  a scaling  humidity  Q* . The 
non-dimensional  parameters  describe  the  fluxes  of  momentum,  heat, 
and  humidity  in  the  surface  layer. 

2.5.1  Moni n-Obukho V Length,  L 


In  a dry  atmosphere  the  characteristic  length 
the  Moni n-Obukho V length. 


DRY 


1 s 


taken  as 


^DRY 


T U. 


g kT, 


(2.6) 


where  T is  mean  surface  layer  temperature,  g is  acceleration  of 
gravity,  and  k is  the  von  Karman  constant.  In  unstable  air 
buoyancy  generated  turbulence  dominates  at  levels  above  -L,  while 
mechanical  turbulence  dominates  much  below  -L.  The  length  is 
corrected  for  moisture  by  using  the  Bowen  ratio  B, 


L 


'-DRY  ^ B^^ 


-1 


(2.7) 


The  Bowen  ratio  is  related  to  the  heat  and  moisture  fluxes  and  is 
expressed  in  terms  of  the  scaling  temperature  and  humidity. 


B 


CnT* 

-P— 

L„Q* 


(2.8) 


where  Cp  is  specific  heat  of  air  at  constant  pressure  and  L^  is 
latent  heat  of  vaporization. 

2.5.2  Scaling  Wind  Speed  U*,  Temperature  T* , and  Humidity  Q* 

The  surface  layer  is  characterized  by  fluxes  of  momentum,  heat 
and  moisture  which  are  independent  of  height.  These  constant 
fluxes  are  used  to  define  scaling  velocity  U* , temperature  T* , and 
humidity  Q* , 


T'w  ' 


U*  = >|-U'w'  (2.9) 

Q*  = - 07  ^ 

where  U',  T'  and  q'  are  fluctuations  of  horizontal  wind,  tempera- 
ture and  humidity  about  their  means,  and  the  overbar  is  an  average 
among  all  fluctuations.  These  scaling  variables  are  related  to 
the  vertical  profiles  of  each  variable, 

U*  = kU[ln(^)  - u.„(z/L)]'’ 

0 

Q*  = |(q-qJCln(^)  - .(/.(z/L ) I"''  (2.10) 

^ ^oQ 

T*  = 

This  is  the  first  time  that  bulk  parameters  are  introduced. 

The  necessary  bulk  parameters  are  wind  speed  U,  normally  taken  at 
10  m height;  air  temperature  T,  again  at  10  m,  air  temperature  at 
the  surface  T^ , taken  to  be  the  sea  surface  temperature;  q, 
specific  humidity  at  10  m;  and  q , specific  humidity  at  the  sea 
surface.  Together  with  these  bulk  parameters,  newly  derived 
parameters  have  been  introduced.  The  and  Z^q  are  rough- 

ness parameters  describing  molecular  transport  of  wind,  temperature, 
and  moisture  near  the  surface,  and  and  are  empirical 

profile  functions.  A new  constant  R is  included  in  the  equations 
for  temperature  and  humidity.  The  value  of  this  constant  is 
accepted  as  0.74  based  on  reliable  field  experiments. 

2.6  PROFILE  FUNCTIONS 

We  adopt  the  profile  functions,  ij; , as  formulated  by  Businger 
et  al  . (1971).  These  functions  have  been  verified  both  over  land 
and  over  water.  These  functions  are 
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L<0 

^ ^ tan(x)  + Tr/2 

X = (1-Y^  z/L)^^ 

H-H  = t^p  = 2 

Constants  in  these  formulations  are  determined  by  fitting  the 
formulations  to  data,  and  correspond  best  to  one  particular  data 
set.  Universal  constants  have  so  far  not  been  found. 

The  hypothesis  that  the  <i/  are  the  same  for  heat  and  moisture 
is  based  mostly  on  measurements  over  land.  There  is  considerable 
controversy  whether  the  assumption  is  valid  over  the  ocean;  but  no 
better  hypotheses  are  available.  Fortunately,  in  all  equations, 
the  Tp's  are  small  compared  to  the  logarithmic  terms,  so  that 
extreme  accuracy  in  the  t/i's  is  not  required.  Further,  ipp  occurs 
only  in  the  correction  of  L for  moisture,  which  itself  is  usually 
a small  correction. 

2.7  ROUGHNESS  LENGTHS  , Z„,.  Z n 

0 0 1 0 U 

The  molecular  transport  of  momentum,  heat  and  moisture  near 
the  surface  is  described  by  the  values  of  Z^ , Z^.|.,  and  Z^p 
respectively.  In  the  marine  environment,  the  roughness  lengths 

p 

have  values  on  the  order  of  1 mm  or  less.  The  dependence  of  C 

7 n 

on  Z is  included  primarily  in  the  function  Iny^,  where  Z is 
elevation,  roughly  10  or  20  meters.  Therefore  is  not  very 
sensitive  to  errors  in  Z^,  although  this  point  should  be  verified. 

The  simplest  assumption  is  that  all  these  roughness  lengths 
are  equal.  It  was  judged  by  the  workshop  participants  that  this 
assumption  is  too  limiting,  and  separate  treatments  of  Z^ , Z^j 
and  Z^p  are  necessary.  One  possible  approach  is  the  technique 
developed  by  Liu  (1978).  Liu  has  derived  relations  between  U and 
U*  and  the  transfer  coefficients  for  momentum,  heat,  and  moisture. 
Transfer  coefficients  for  momentum,  heat  (Cp,  Cp),  and  moisture 
are  directly  related  to  the  roughness  lengths. 


L>0 


(2.11) 


'^H‘’''Q  ■■  R L 
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ln(z/Zj  = ^ + -i.^,(z/L) 


InCz/Z^j)  " R iI^hCz/L) 


(2.12) 


where  z is  10  tn.  The  derived  relationships  are  displayed  in 
Fig.  2.1.  The  values  of  the  roughness  lengths  cannot  be  obtained 
in  closed  form  expressions,  and  the  following  technique  is 
recommended  : 

a.  Z^ , Z^g  at  neutral  are  found  directly  from  the 

curves  of  Fig.  2.1,  using  the  value  of  10  m wind. 

b.  After  stability  corrections  and  are  calculated,  a new 
set  of  roughness  lengths  is  found  from  Figure  2.1  using  the 
roughness  Reynolds  number. 


R 


r 


(2.13) 


where  v is  the  kinematic  viscosity  of  air. 

2.8  SEA  SURFACE  TEMPERATURE,  T^ 

2 

The  Cj  can  now  be  evaluated  in  terms  of  bulk  meteorological 
parameters,  and  one  other  parameter  surface  temperature,  T^ , not 
normal  1 y observed  . 

In  most  cases,  bucket  temperature,  T^,  rather  than  the  surface 
temperature  is  Pleasured.  We  estimate  T^-T^  from  the  theory  by  Liu 
( 1 978).  Liu,  (Fig.  2.2),  gives  T -T  as  a function  of  the  wind 

^ W 

speed  and  the  difference  between  bucket  temperature  and  air 
temperature  at  10  m height,  aT.  In  practice,  temperatures  at  any 
height  between  7 m and  20  m should  be  satisfactory  in  place  of  the 
temperature  at  10  m. 

Liu's  diagram  is  approximated  by  a set  of  equations  which 
should  be  more  usable  for  computation  than  his  figure.  These 
equations  are: 
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2.1.  Bulk  Transfer  Coefficients  at  Neutral  Conditions  Versus 
10  m Wind  U and  Roughness  Reynolds  Number  R^  (Liu,  1978). 

The  transfer  coefficients  are  approximately  calculated  from 
wind  U;  when  scaling  speed  U*  and  roughness  length  Z are 
found,  transfer  coefficients  are  found  from  R usi ng°formulas 
in  text. 


Fig.  2.2  Temperature  Difference  Across  the  Aqueous  Sublayer  Versus 
10  m Wind  for  Values  of  Air-Bulk  Water  Temperature  Difference 
Ty^-Tio  (Liu,  1 978).  These  formulas  account  for  radiation  and 
molecular  effects  in  the  surface  layer. 


(2.14) 


2.9  EVALUATION  OF  CONSTANTS 

Certain  empirical  constants  were  introduced  in  the  definition 

2 

of  scaled  quantities,  the  vertical  flux  profiles,  and  in  the  Cy 

stability  function.  These  constants  depend  on  the  particular  data 

set  from  which  they  were  determined,  and  there  is  no  uniform 

agreement  in  the  literature  on  the  best  values  of  each.  In 

choosing  the  constants,  a set  of  values  corresponding  to  one 

particular  data  set  should  be  chosen  because  the  constants  are 

interdependent.  Table  1 provides  the  ranges  of  the  various 

2 

constants  for  use  in  a C^  model  together  with  the  recommended 
value  that  emerged  from  the  NEPRF  workshop. 


Constant 

Table  1 

Range 

Recommendati on 

For 

Scaling  Functions 

k 

0 . 33-0.41 

0.37 

R 

.74-1 

.74 

For 

Profile  Functions 

^m 

15-16 

1 5 

9-18 

9 

4.7 

4.7 

For 

Stability  Function 

b 

.667-1 

1 . 

2.10  CALCULATION  OF  C^ 

n 

2 

The  calculation  of  C^  proceeds  using  the  equations  described 
in  Paras.  2.2  through  2.8  and  the  constants  given  in  Para.  2.9. 
The  step-by-step  procedure  is  as  follows: 
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a.  Modify  the  sea  surface  temperature  measurement  to  an 
actual  sea  surface  temperature  value  using  Eq.  (2.14)  and  T,  T-T^ 
and  U . 

b.  Set  Initial  stability  to  neutral,  i.e.,  z/L  = 0 and 

c.  The  bulk  wind  speed  determines  an  initial  value  of 

Cp  . , Cq  from  Fi  g . 2.1. 

d.  Calculate  , Z^j,  and  Z^p  from  Cp,  Cj^,  Cp  with 
Eq.  (2.12). 

e.  Evaluate  scaling  parameters  U*,  T*  and  Q*  using 
Eq.  (2.10). 

f.  Evaluate  Bowen  ratio,  B,  using  Eq.  (2.8). 

g.  Evaluate  Mon i n-Obukho v length  for  moist  atmosphere  using 
Eqs . (2.6)  and  (2.7). 

h.  Test  if  the  new  L is  within  1%  of  old  value  of  L.  If  so, 
go  to  statement  1.  If  not,  go  to  statement  i. 

i.  Calculate  stability  functions  using  Eq.  (2.11). 

j.  Calculate  rougnness  Reynolds  number  Rr  from  Eq.  (2.13). 

k.  Evaluate  Cp,  C^,  Cp  from  Fig.  2.1  using  Rr  and  go  to 
statement  d. 

2 

l.  Evaluate  Cj  stability  function  (g(z/L)  using  Eq.  (2.5). 

2 

m.  Evaluate  C,  using  Eq.  (2.4). 

2 

n.  Calculate  C^  using  Eq.  (2.2). 

This  procedure  is  graphically  outlined  in  Fig.  2.3. 

2.11  DEVIATION  OF  C^  FROM  ENSEMBLE  MEAN  VALUE 

2 

The  values  of  C^  predicted  by  the  surface  model  are  values  of 
mean  turbulent  quantities.  The  equations  were  derived  from 
measurements  of  temperature  fluctuations  of  one  turbulent  eddy 
size  at  one  point  averaged  over  long  periods  of  time,  i.e.,  20  to 
30  minutes.  The  laser  system,  on  the  other  hand,  is  affected  by 
turbulence  over  a spectrum  of  eddy  sizes  and  operates  over  a short 
time  period,  about  10  seconds,  and  over  a long  path,  about  10  km. 
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The  averaging  time  to  obtain  mean  fluctuations  of  refractive 
index  is  proportional  to  the  relevant  eddy  size,  and  inversely 
proportional  to  the  wind  speed.  In  the  case  of  the  measurements 
on  which  the  scaling  laws  are  based,  the  eddy  size  corresponds  to 
a temperature  probe  separation  r,  about  20  cm.  Therefore  the 
averaging  time,  T^ , relevant  to  the  scaling  laws  is 

Tg  = r/U  . (2.15) 

In  the  case  of  an  optical  system,  the  eddy  size  corresponds 
to  a combination  of  wavelength  \ and  range  R.  Therefore  averaging 
time  T^p  for  the  optical  system  is 

Tqp  = /WU  • (2.16) 

2 

Wyngaard  and  Clifford  (1978)  show  that  the  mean  corres- 
ponding to  the  two  methods  are  the  same  if 

2 

From  this  we  infer  that  a laser  beam  is  affected  by  a mean 

corresponding  to  a model  prediction  if  the  minimum  time  interval 

of  the  laser  transmission  T- ^ satisfies 

T r 

T.^^  = /XR/144  • (2.18) 

For  a laser  of  wavelength  3.5  ym,  and  range  10  km,  and  a model 
ensemble  averaging  time  of  20  minutes,  the  equation  states  that 
the  turbulence  model  prediction  of  describes  the  turbulence 
affecting  the  system  if  the  laser  transmission  duration  is  at 
least  1.5  sec  . 

2.12  RECOMMENDATIONS  - SURFACE  LAYER  C^  ASSESSMENT 

2 

The  surface  layer  C model  outlined  here  is  recommended  for 

n 2 

verification  and  use  by  the  HEL  program  to  predict  C^  for  within 
the  atmospheric  surface  layer  (within  100  m of  the  sea  surface). 
The  model  is  basically  the  NEPRF  preliminary  model  (Appendix  C) 
with  the  following  modifications: 


a.  Model  should  account  for  humidity  using  Eq.  {2.2). 

b.  Roughness  lengths  for  momentum,  heat,  and  moisture 
should  be  treated  as  different  functions,  as  discussed  in 
Para  . 2.7  . 

c.  Sea  surface  temperature  should  be  derived  from  bucket 
sea  surface  temperature  measurement  using  Eq.  (2.14). 

The  humidity  modification  has  not  been  verified  in  the  marine 

atmosphere.  A limited  field  program  should  be  conducted  to  verify 

2 2 

the  recommended  modification  of  Cy  to  produce  C^. 
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3.  REFRACTIVE  INDEX  STRUCTURE  PARAMETERS:  PROFILES  IN  THE 
UPPER  MARINE  PLANETARY  BOUNDARY  LAYER 


3.1  INTRODUCTION 

We  now  discuss  methods  of  computing  vertical  profiles  of 
refractive  index  structure  parameters  above  the  surface  layer. 

In  particular,  this  chapter  is  concerned  with  that  portion  of  the 
atmosphere  between  the  top  of  the  surface  layer  and  the  top  of  the 
marine  planetary  boundary  layer  (MPBL).  We  term  this  atmospheric 
region  the  upper  MPBL. 

Chapter  2 presents  well-established  Moni n-Obukhov  scaling 
2 

laws  that  yield  C^  within  the  surface  boundary  layer.  The  scaling 
laws  change,  however,  as  one  moves  to  higher  altitudes.  The 
mechanical  surface  turbulence  becomes  less  important  than  surface 
heating  and  the  influence  of  the  inversion.  The  stability  of  the 
MPBL,  furthermore,  may  profoundly  alter  the  nature  of  the  scaling 
laws  in  the  upper  MPBL. 

When  the  MPBL  is  stable  (L>0),  no  adequate  theory  exists  for 

scaling  atmospheric  variables  in  the  upper  MPBL.  This  is  an  area 

of  active  research;  new  theories  can  be  expected  to  emerge.  When 

the  MPBL  is  unstable  (L<0),  there  are  scaling  laws  for  the  upper 

MPBL.  These  scaling  laws,  however,  require  knowledge  of  universal 

functions  analogous  to  the  universal  function  g(z/L)  that  occurs 

2 

in  the  expression  for  Cj  in  the  surface  layer.  These  universal 

functions  in  the  upper  MPBL  are  unknown  currently.  The  scaling 

laws  for  the  upper  MPBL  remain  qualitative. 

Upper  MPBL  scaling  is  detailed  in  Para.  3.2.  In  Para.  3.3  we 

discuss  numerical  MPBL  models  and  how  they  can  contribute  to  the 

study  of  structure  parameter  profiles.  Observational  data, 

primarily  from  instrumented  aircraft,  are  discussed  in  the 

following  section.  In  the  final  section,  we  recommend  methods  of 

2 

obtaining  a better  understanding  of  in  the  upper  MPBL. 
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3.2  SCALING  LAWS  IN  THE  UPPER  MPBL  UNDER  UNSTABLE  CONDITIONS 
3.2.1  Regions  of  Val 1 di ty 


The  depth  of  the  MPBL  Is  generally  defined  by  the  Inversion 
height,  Z. . In  the  surface  layer,  where  Z<0.1Z^,  Moni n-Obukho v 
scaling  has  been  adequately  confirmed  under  a variety  of  circum- 
stances. Near  the  Inversion,  mixed  layer  scaling  Is  appropriate. 
In  this  region,  where  Z.>Z>-L,  the  velocity  scale  analogous  to  U* 
In  the  surface  layer  Is  W*  and  the  temperature  scale  analogous  to 
T*  In  the  surface  layer  Is  9*,  where 


and  H Is  the  surface  heat  flux;  p,  atmospheric  density;  and  Cp, 
the  specific  heat  at  constant  pressure. 

Between  the  two  regions,  the  nature  of  the  scaling  depends  on 
the  strength  of  shear  production  relative  to  the  strength  of 
buoyant  production  of  turbulence.  When  buoyant  production  Is 
dominant  In  the  surface  layer,  characterized  by  -L<0.1Z^.  , the  two 
scaling  regions  overlap.  The  overlap  region,  -L<Z<.1Z^.  , Is  termed 
a layer  of  free  convection.  Here,  Moni n-Obukhov  and  mixed  layer 
scaling  apply  simultaneously. 

In  the  other  case,  where  shear  production  of  turbulence  Is 
not  negligible  compared  to  buoyant  production  within  the  surface 
layer,  I.e.,  -L>0.1Z^,  no  overlap  region  exists.  Rather,  the 
layer  0.1Z^<Z<-L  Is  a transition  region  between  the  Moni n-Obukho v 
and  mixed  layer  scaling  regions  where  neither  scaling  law  applies. 
This  region  Is  an  enigma  Insofar  as  scaling  Is  concerned.  A MPBL 
having  a free  convective  region  and  a MPBL  having  a transition 
layer  are  shown  In  the  Fig.  3.1. 
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(a) 


Mixed  Layer  Scaling 


M-0  J } Free  Convection 
Scaling 

4 


Mixed  Layer  Scaling 

Transition  Region 
M-0  Scaling 


Fig.  3.1.  (a)  Boundary  layer  which  has  a free  convection  region 
in  which  both  mixed  layer  scaling  and  Moni n-Obukhov  (M-0) 
scaling  apply  simultaneously;  and  (b)  a boundary  layer  with 
no  free  convection  region.  Neither  mixed  layer  scaling  nor 
M-0  scaling  apply  in  the  transition  region. 
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3.2.2  Mixed  Layer 

In  the  region  Z^.>Z>-L,  mixed  layer  scaling  requires  that 


, (3.2) 

where  F(Z/Z^) is  a function  that  has  not  been  determined.  Note 
that  in  mixed  layer  scaling,  e*  and  Z^  are  the  relevant  scaling 
parameters.  Frisch  and  Ochs  (1  975  ) attempt  to  evaluate  F(Z/Z^.) 
by  using  aircraft  data  collected  in  the  marine  environment  near 
San  Diego.  They  find 


F(Z/Z^)  = 1 .0+0.84Z/Z^+4.13(Z/Z.)^  , 
0<Z/Z^<0.8  . 


(3.3) 


The  data  upon  which  this  expression  is  based  contains  considerable 
scatter,  particularly  in  the  upper  half  of  the  MPBL  . Equation 
(3.2)  and  (3.3)  must  be  applied  with  caution.  In  the  mixed  layer, 
the  properties  of  the  free  atmosphere  above  the  mixed  layer  may 
also  become  important.  For  instance,  Deardorff  (1974)  finds  that 
entrainment  at  the  top  of  a growing  mixed  layer  can  significantly 
influence  turbulence  statistics  involving  humidity  throughout 
much  of  the  boundary  layer.  Considerable  generalization  of 
Eq.  (3.2)  may  be  required  before  these  types  of  influences  can  be 
handled  by  scaling  laws. 

3.2.3  Free  Convection  Scaling 

2 

Within  a layer  of  free  convection  -L<Z<0.1Z.,  C,  is  not 

1 1 2 

influenced  by  shear  turbulence  production,  that  is,  Cj  is 

independent  of  U* . This  condition  is  satisfied  if  g(z/L)  is 

-2/3 

proportional  to  Z ' , in  which  case 

C^  * AZ''^^^H^^^(g/T)"^^^  , (3.4) 

where  A is  a constant.  Within  the  region  of  applicability,  this 
relation  fits  observations. 
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3.2.4  Transition  Region 

When  shear  production  is  not  negligible  anywhere  within  the 
surface  layer,  -L>0.1Z.,  scaling  is  complicated.  The  middle 
region  of  the  MPBL  is  influenced  both  by  surface  shear  turbulence 
production  and  the  presence  of  the  inversion.  Here  U*  and  Z.  are 
important  scaling  parameters,  and  we  must  write 


= Z/L)  , 


(3.5) 


where  F'  is  still  another  universal  function.  This  scaling  law 
transition  region  is  quite  common  over  the  oceans.  No  quantita- 
tive proposal  for  the  form  of  F'  exists. 


3.3  NUMERICAL  MPBL  MODELS:  CALCULATION  OF  STRUCTURE  PARAMETERS 
3.3.1  Basic  Nature  of  Models 

Numerical  boundary  layer  modeling  involves  solving  equations 
which  describe  the  basic  laws  of  fluid  dynamics  and  atmospheric 
thermodynamics.  These  equations  generally  can  be  characterized 
by  the  expression 

•||-  = Source/sink  + mean  transport  + turbulent  transport 

where  41  is  any  one  of  the  dependent  variables.  For  example,  if  41 
represents  temperature,  then  the  source/s i nk. term  would  involve 
processes  such  as  radiative  transfer,  latent  heat  release,  etc. 

The  mean  transport  term  would  represent  temperature  advection  by 
the  mean  wind,  and  turbulence  transport  would  involve  the  diver- 
gence of  the  turbulent  heat  flux. 

The  source/sink  terms  and  the  mean  transport  terms  in  the 
fundamental  equations  are  generally  treated  similarly  in  all 
models.  Substantial  differences  exist,  however,  in  the  manner  in 
which  models  treat  the  turbulent  transport  terms.  We  will  discuss 
these  different  treatments  of  the  turbulent  transport  terms  in 
more  detail  in  this  section.  Here  we  note  only  that  since  our 
concern  is  with  turbulent  refractive  index  fluctuations,  substan- 
tial attention  and  care  must  be  given  to  the  handling  of  the 
turbulent  transport  terms. 
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Our  interest  here  centers  upon  two  different  types  of 
numerical  MPBL  models.  One  type  we  term  a subgrid-scale  closure 
model  while  the  other  type  we  term  a second-moment  turbulent 
cl osure  model  .* 

3.3.2  Subgrid-Scale  Closure  Models 

Subgrid-scale  closure  models  use  the  fundamental  equations 
to  describe  all  the  turbulent  eddies  within  the  resolution  of  the 
grid.  The  smaller  subgrid-scale  eddies  must  be  parameterized. 

An  example  of  a subgrid-scale  closure  model  is  that  developed 
by  Deardorff  (1973).  In  his  model,  the  parameterized  subgrid- 
scale  turbulent  fluxes  generally  represent  small  fractions  of  the 
total  fluxes.  Thus,  the  model  forecasts  are  not  very  sensitive 
to  the  manner  in  which  the  subgrid-scale  is  parameterized. 

Because  of  these  considerations,  the  subgrid-scale  closure 
model  has  considerable  appeal  as  a tool  for  investigating  basic 
physical  processes  in  the  MPBL.  The  major  disadvantage  inherent 
in  this  type  of  model  involves  the  computational  requirements. 

For  example,  a version  of  the  Deardorff  model  used  to  simulate  an 
undisturbed  trade  wind  MPBL  took  about  30  h computer  time  with  a 
CDC  7600  to  produce  a 5 h model  time  simulation  (Sommeria,  1976). 
This  does  not  include  any  development,  debugging,  or  preliminary 
testing  time.  The  time  to  analyze  the  model  results,  furthermore, 
is  comparable  to  that  required  in  an  extensive  field  experiment. 
Therefore,  it  seems  prohibitive  at  this  time  to  examine  in  a 
wide  variety  of  atmospheric  situations  with  a subgrid-scale 
closure  model  . 


*The  term  'closure'  refers  to  the  process  whereby  a set  of 
equations  becomes  closed,  that  is,  the  set  has  an  equal  number  of 
equations  and  unknowns.  We  find  "second-moment"  closure  to  be 
more  descriptive  than  "second-order"  or  "higher-order"  closure 
because  the  assumptions  that  close  the  set  are  introduced  into  the 
equations  for  the  second  moments. 
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3.3.3  Second-Moment  Closure  Models 

Second-moment  closure  models  solve  the  fundamental  equations 
for  the  total  turbulent  flux  rather  than  resolving  a portion  of 
this  flux  by  explicitly  computing  the  turbulent  flow  on  scales 
larger  than  the  grid  scale.  The  parameterization  involved  in  the 
closure  assumptions  is  more  important  here  than  in  subgrid-scale 
closure  models;  the  parameterization  affects  the  total  flux,  not 
just  that  portion  of  the  flux  that  is  below  the  grid  resolution. 

The  treatment  of  the  total  turbulent  flux  can  involve  parameter- 
ization at  several  levels  corresponding  to  the  level  of  the 
moment  closure  assumptions. 

These  parameteri zations  are  best  seen  by  first  considering 
first-moment  closure,  or  K-theory.  Here  the  fluxes  of  momentum, 
heat,  and  moisture  that  appear  in  the  fundamental  equations  are 
replaced,  in  analogy  with  molecular  diffusion,  with  terms  involving 
the  product  of  an  eddy  coefficient  and  a mean  gradient  of  wind, 
temperature,  or  humidity.  This  is  first-moment  closure  because 
the  assumptions  which  yield  a closed  set  are  introduced  into  the 
equations  for  the  first  moments,  that  is,  the  means.  First-moment 
closure  has  serious  shortcomings,  however,  because  in  actuality 
fluxes  are  not  so  clo,sely  tied  to  the  local  gradients.  The  eddy 
coefficient,  furthermore,  is  a function  of  the  flow  itself  and, 
therefore,  cannot  a priori  be  specified. 

In  second-moment  closure,  the  turbulent  fluxes  (the  second 
moments  of  the  flow)  in  the  fundamental  equations  are  not  directly 
replaced  as  in  first-moment  closure.  Rather,  time-dependent 
equations  describing  the  dynamic  behavior  of  these  fluxes  are 
developed.  These  equations  contain  additional  higher-order  moment 
terms.  These  unknowns  are  parameterized  in  terms  of  the  known 
variables,  thereby  closing  the  set  of  equations.  Thus,  in  a 
second-moment  closure  model,  one  charts  the  progress  of  both  the 
mean  variables  and  the  turbulence  moments. 

The  upper  boundary  conditions  of  these  numerical  models  are 
usually  selected  to  be  in  the  free  atmosphere  where  turbulent 
fluxes  are  weak.  Near  the  surface,  Moni n-Obukho v similarity  theory 
provides  boundary  values.  Thus,  a surface  layer  model  as  discussed 
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in  Chapter  2 Is  Incorporated  as  a part  of  the  lower  boundary 
condition  in  these  numerical  models. 


2 

3.3.4  Calculating  C^  with  a Second-Moment  Closure  Model 

Although  sensitive  to  the  closure  assumptions,  second-moment 

closure  models  can  answer  some  of  the  questions  concerning  the 
2 

nature  of  C^  profiles  throughout  the  MPBL  . As  mentioned  in  Chapter 
2,  Wesely  (1976)  shows  that  the  optical  structure  parameter  may  be 
written  as 


C 


2 
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2(1-A2/A^ )T 
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(l-A^/A, )T 
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If  this  equation  is  to  be  useful,  we  need  a way  of  computing  Ct, 

2 * 
Cgj,  and  Cg  with  the  closure  model  throughout  the  MPBL.  We  now 

discuss  a possible  approach. 

In  terms  of  dissipation  rates,  the  temperature  and  humidity 
structure  parameters  may  be  written  (Wyngaard  et  al . , 1978) 


= ABEgC 


-1/3, 


'eT 


-1/3 


(3.7) 


Here  e is  the  turbulent  kinetic  energy  dissipation  rate;  zj  and 
Eg  are  the  dissipation  rates  of  temperature  and  vapor  pressure 
variances,  respectively;  Egj  is  the  dissipation  rate  of  temperature- 
humidity  covariance;  and  b and  y are  constants. 

The  terms  e-j-.  Eg,  Egj  and  e appearing  in  these  structure 
parameter  relationships,  also  appear  in  the  equations  that  form 
the  core  of  a second-moment  closure  model.  These  dissipation  terms 
are  approximated  by  closure  assumptions.  The  typical  closure 
assumptions  for  dissipation  terms  are 


E 


^T  “ 


(3.8) 


where  T = )i/q  is  a turbulent  time  scale  and  t is  a turbulent  length 
scale.  Here  is  twice  the  turbulent  kinetic  energy  (i.e., 

2 2 2 2 

(i.e.,  q = u'  +v ' +w ' ).  Using  these  expressions  the  structure 
parameters  are 


.2/3  ’ 


e '2 

“^2/3  ’ 


(3.9) 


Second-moment  closure  models  compute  the  evolution  of  T'^,  e'^, 

2 2 

T'e',  and  i.  Thus  we  solve  for  C^,  C^,  and  using  Eq . (3.9) 

and  for  optical  using  Eq.  (3.6). 

n 2 

Computation  of  throughout  the  MPBL  as  we  have  been 

describing  is  the  subject  of  a recent  report  by  Burk  (1978).  He 

uses  a one-dimensional  second-moment  closure  model.  The  question 

of  what  dimensionality  a model  should  have  to  adequately  describe 
2 

in  the  MPBL  is  the  subject  of  the  next  subsection. 

3.3.5  Multidimension  vs.  1-D  Models 

The  dimensionality  of  a second-moment  closure  model  can  also 
be  an  important  consideration  in  seeking  vertical  structure 
parameter  profiles.  A 1-D  model  assumes  horizontal  homogenity. 

It  cannot  yield  information  concerning  the  nature  of  large  two  or 
three  dimensional  turbulent  structures  such  as  horizontal  roll 
vortices  or  hexagonal  cells.  For  this  type  of  information,  a 
multidimensional  model  is  required. 

Teske  and  Lewellen  (1979)  use  a two-dimensional  second-moment 
closure  model  to  study  horizontal  roll  vortices.  They  find  that 
the  upward  moving  region  of  a roll  concentrates  the  turbulent 
fluctuations.  This  is  particularly  evident  near  the  top  of  the 
boundary  layer.  The  two-dimensional  model  also  can  compute  wave 
motion  near  an  inversion,  which  a 1-D  model  cannot.  As  a result, 
Teske  and  Lewellen  found  the  temperature  variance  near  the  inversion 
computed  with  the  two-dimensional  model  to  be  about  twice  that  found 
with  a 1-D  model.  This  may  be  a unique  result  and  should  not  be 
used  for  generalizing  about  the  advantages  of  dimension. 


One  further  practical  consideration  enters  the  picture  when 
discussing  the  utility  of  1-D  vs.  multidimensional  models.  In  the 
marine  environment,  one  rarely  has  the  type  of  high  resolution 
data  necessary  to  provide  initial  and  boundary  conditions  to  a 
multidimensional  MPBL  model.  Typically,  only  isolated  radiosonde 
and  bulk  air-sea  ship  data  area  available.  Such  a sparse  data  base 
is  more  compatible  with  the  initial  and  boundary  condition  require- 
ments of  a 1-D  model.  Only  further  study  will  provide  answers 
concerning  the  relative  merits  of  1-D  vs.  multidimensional  models. 

3.4  OBSERVATIONS  AND  VERIFICATIONS 

A numerical  model  can  aid  in  our  understanding  of  structure 
parameter  profiles.  But  observational  data  must  be  continually 
used  as  a touchstone.  Idealized  model  computations  can  have  an 
elegant  simplicity  and  beauty  that  fades  when  real  data  problems 
must  be  addressed.  Observations  provide  the  only  means  whereby  a 
numerical  model  can  be  improved  and  validated. 

Measurements  of  structure  parameters  in  the  marine  environment 

are  limited.  The  Cooperative  Experiment  in  West  Coast  Oceanography 

2 

and  Meteorology  (CEWCOM)  contains  a set  of  aircraft  C^  measurements 
in  conjunction  with  bulk  air-sea  values.  Such  measurement  programs 
remain  an  active  area  of  research  at  the  Naval  Postgraduate  School 
(Fairall  et  al.,  1978;  Davidson  et  al.,  1978).  Preliminary  results 
from  a recent  field  experiment  off  the  coast  of  Panama  City, 

Florida,  were  reported  on  at  this  workshop  by  Dr.  Fairall.  This 
data,  and  that  from  other  large  research  experiments  such  as  AMTEX 
and  GATE,  need  to  be  exploited  more  fully  both  for  model  validation 
purposes  and  for  the  purpose  of  testing  scaling  laws. 

Care  must  be  exercised  when  aircraft  structure  parameter 
measurements  are  compared  with  closure  model  results  or  with 
scaling  law  predictions.  Scaling  laws  and  closure  models  deal,  in 
principle,  with  ensemble-averaged  turbulence  statistics.  As  we 
know,  however,  ensemble  averaging  is  hardly  practical  in  geophysical 
flows,  and  instead  one  generally  averages  atmospheric  measurements 
over  time.  The  prime  question  then  is:  what  ave  aging  time  of 
aircraft  data  is  necessary  to  obtain  statistics  which  are  relevant 
to  the  predictions  of  the  models  or  the  scaling  laws? 


The  standard  two-pronged  temperature  probe  measurements  of 

2 

generally  require  about  30  minute  averages  to  obtain  stable 
statistics  when  the  measurements  are  made  from  a stationary  plat- 
form, If  the  instrument  is  attached  to  an  aircraft,  about  an 
order  of  magnitude  advantage  in  averaging  time  can  be  obtained. 
Thus,  a 3 minute  average  at  each  flight  level  corresponds  in 
accuracy  to  the  30  minute  average  from  a stationary  platform. 
Shorter  time  averaging  of  the  aircraft  data  will  yield  profiles 
having  a very  jagged  appearance.  These  short-time  averaged 
profiles  cannot  be  directly  related  to  scaling  law  or  closure 
model  results. 

3.5  RECOMMENDATIONS  - UPPER  MPBL  ASSESSMENT 

a . Verify  1-D  Closure  Model  with  Existing  Data 

2 

The  first  step  towards  understanding  C^  in  the  upper  MPBL  by 
means  of  a numerical  model  must  involve  real  data  verification  of 
the  model.  The  existing  volume  of  data  must  be  searched  for  those 
cases  in  which  sufficient  information  is  available  to  supply  the 
initial  and  boundary  condition  requirements  of  the  model.  Then 
comparison  tests  of  model  predictions  and  measured  structure 
parameter  profiles  can  be  made.  Preliminary  work  along  these  lines 
is  reported  by  Burk  (1978).  These  efforts  are  continuing. 

Primarily  for  reasons  of  computational  economy,  a 1-D  second- 
moment  closure  model  is  recommended  for  this  numerical  model  study. 
When  possible,  comparison  between  1-D  and  multidimensional  model 
results  is  recommended.  Comparisons  among  different  versions  of 
1-D  closure  models  are  also  desirable.  Particular  attention  should 
be  paid  to  improving,  where  possible,  the  closure  assumptions  which 
are  used  in  these  models. 

b . Use  Closure  Models  to  Derive  Scaling  Laws 

As  noted  previously,  if  scaling  laws  are  to  be  useful,  the 
functional  form  of  the  universal  functions  must  be  known.  For 
example,  in  Eq.  (3.2)  we  need  to  know  F and  in  Eq.  (3.5)  we  need  to 
know  F'.  We  recognize  that  these  functions  may  depend  on  the 
entrainment  of  temperature  and  moisture  from  the  atmosphere  above 
the  inversion. 


35 


The  universal  functions  can  be  established  either  from 
observations  or  from  mathematical  modeling.  The  latter  technique 
is  much  more  flexible  in  that  many  alternatives  can  be  studied 
with  a single  program.  Hence,  we  recommend  that  a numerical  model 
be  used  to  help  develop  general  scaling  laws  which  describe  in 
the  upper  MPBL. 


c . Conduct  Experimental  Verification 

The  existing  data  base  of  structure  function  measurements  in 
the  MPBL  is  small.  This  is  particularly  true  when  those  cases  are 
discarded  that  have  insufficient  data  for  initializing  a numerical 
model.  For  example,  many  of  the  aircraft  measurements  of  have 
no  wind  measurements,  making  estimates  of  turbulent  shear  produc- 
tion s pecul ati ve  . 

The  existing  data  base  is  inadequate  for  testing  the  predic- 
tions of  a closure  model.  Further,  when  constructing  generalized 
scaling  laws  of  the  upper  MPBL,  new  measurements  are  needed. 
Therefore,  a new  detailed  field  experiment  is  recommended.  This 
field  experiment  is  to  be  limited  in  scope.  The  preliminary 
scaling  law  and  numerical  model  studies  will  define  and  dictate 
the  scope  of  the  new  field  experiment. 
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APPENDIX  C 


PRELIMINARY  TURBULENCE  MODEL 

The  main  atmospheric  effects  on  an  optical  signal  are 
absorption  and  scattering  by  molecules  and  aerosols  and  degrada- 
tion of  the  signal  by  atmospheric  turbulence.  The  following 
presents  a preliminary  model  of  the  intensity  of  atmospheric 
turbulence  as  a function  of  bulk  meteorological  variables.  This 
model  is  currently  used  by  the  Navy  HEL  program.  At  this  point, 
however,  the  model  is  to  serve  as  a start  of  discussion  with  all 
segments  open  to  modification. 

The  effect  of  refractive  index  inhomogeneities  on  a coherent 
laser  beam  is  to  focus  a segment  of  the  beam  to  a point  determined 
by  the  size  and  location  of  the  inhomogeneity.  Over  a long  path 
this  causes  a defocusing  and  loss  of  coherence  in  the  beam.  An 
analysis  of  turbulent  effects  on  beam  characteristics  requires 
assessment  of  a mean  value  of  the  intensity  of  turbulence  over 
the  beam  path,  and  for  the  time  of  transmission. 


Calculation  of  Optical  Parameters 

The  inhomogeneities  of  index  of  refraction  are  described  by 

2 

an  index  of  refraction  structure  constant  C^,  defined  as 


C 


2 

n 


[n(r^ )-n(r2)]^ 


( ki-''2l  ^ 


-2/3 


(1) 


The  optical  index  of  refraction  of  air  is  related  to  tempera 
ture  T and  water  vapor  pressure  e. 


n = 1 + A-j  ^2  T 

where  A^  and  A2  are  constants. 

The  fluctuation  of  refractive  index  are  therefore  related  to 
fluctuations  of  temperature  and  humidity  and  their  correlation 
(Wesely,  1973), 


■ '^eT  * “37 


(3) 


where  B^  , B^.  and  B^  are  constants. 

The  first  term  Is  dominant  in  dry  conditions.  The  second 
term  can  contribute  significantly  under  special  humidity  condi- 
tions described  by  Friehe  (1977).  The  preliminary  model  neglects 
the  latter  two  contributions  for  two  reasons.  First,  the 
contribution  is  the  dominant  conttribution  to  the  refractive  index 

structure  function.  Second,  there  are  no  preliminary  parameteri- 
2 2 

zations  of  and  C,  available  which  have  been  verified  as 
eT  T - 

extensively  as  those  of  (Wesely,  1973;  Friehe,  1977). 

By  a combination  of  dimensional  arguments  and  curve  fitting 

to  the  Kansas  data  base,  Wyngaard  et  al  . (1971)  derived  a semi- 

2 

empirical  relationship  of  C.|.  and  turbulence  parameters. 


(4) 


where 


g(z^/L) 


4. 9(1-7  Z 

d 

4.9(1+2.75  Z /L) 

a 


Za/L  < 0 
Za/L  > 0 


This  parameterization  seems  to  also  be  valid  over  the  ocean,  if 
proper  experimental  precautions  are  taken  during  the  measurement 
(Friehe,  1977). 


Calculation  of  Turbulence  Parameters 

The  turbulence  parameters  T*  and  L are  calculated  from  the 
bulk  meteorological  variables  using  semi -empi ri cal  heat  and 
momentum  flux  profiles.  Two  commonly  used  formulations  of  these 
profiles  are  the  KEYPS  (Lumley  and  Panofsky,  1964)  and  the 
Businger-Dyer  formulas  (Businger  et  al  . , 1971).  The  Busi nger-Dyer 
formalism  was  chosen  to  represent  the  flux  profiles  because  it  is 
widely  used  and  is  part  of  the  approximation  formulas  to  be 
described  below. 
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The  calculation  of  the  height  by  the  Moni n-Obukho v length 
Z^/L  proceeds  with  the  approximation  used  by  Barker  and  Baxter 

d 

(1975).  The  scaled  height  is  given  by 


U (T  +T  ) 
a'  va  VO ' 


where  (for  unstable  conditions) 


'*'m  " ■ )(^)  ] + 2arctan(e)  j 


c = [1-v,  Z^/L]-- 

-Ph  = - 2 In  [^(1  + /1-Yh  Zg/L)] 

Using  the  definition  of  bulk  Richardson  number 

2gZ  T -T 

R = a va  VO 

•^i  „ 2 T +T 

V.  va  VO 

d 

z 

Barker  and  Baxter  find  an  approximate  form  for  , 

^ = Ri-[0.47li  In(^)  - 1 .045]  (7) 

0 

This  equation  is  inaccurate  for  0<^£-0.05  in  which  case  one  uses 


«i  - I? 


Mr  + ^ 
e 1 
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1-6R. 


The  temperature  function  T*  is  obtained  from  an  integration 
of  the  flux  profile  relationships 
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(9) 


■ 


T. 


R[1  n(Y^)+t()^ 


The  several  constants  in  ‘ihe  flux  profile  relationships  have 

a range  of  values  attributed  to  them  in  the  literature.  Yaglom 

(1974)  points  out  that  the  momentum  profile  constant  v is 

m 

accepted  at  about  15  or  16  by  several  authors.  The  heat  profile 
constant,  however,  has  been  given  a range  of  values  from  9 
(Businger  et  al.,  1971)  from  the  Kansas  data  to  18  over  the  Utah 
desert  (Mon ji  , 1 973). 


Laboratory  measurements  of  the  von  Karmann  constant  k range 
from  0.2  to  0.8,  although  atmospheric  measurements  show  values 
from  0.3  to  0.44.  The  roughness  length  over  the  ocean  is  generally 
determined  from  Charnock's  (1955)  relation 


(10) 


where  ot  ranges  from  0.077  (Roll,  1 965)  to  0.0144  (Garrett,  1 977  ). 

In  choosing  the  values  of  the  constants  for  this  model,  it 
was  realized  that  there  is  an  interdependence  among  them  during 
most  measurement  programs.  Hence  we  chose  a set  of  constants 
consistent  with  one  experimental  set  of  data.  Since  the  Kansas 
data  is  widely  accepted  as  a careful  and  extensive  set  of  measure- 
ments, the  constants  developed  by  Businger  et  al  . (1971)  are  being 
used.  These  are 


k = 0.35,  Yjjj  = 16,  = 9 


Because  of  the  uncertainty  in  the  calculation  of  Z , a 

-4  ° 

constant  value  of  4x10  m was  chosen,  independent  of  wind  speed. 


GLOSSARY  FOR  PRELIMINARY  TURBULENCE  MODEL 

2 

Cj  - Temperature  structure  function 
2 

- Index  of  refraction  structure  function 
n 

2 

Cg  - Water  vapor  density  structure  function 
2 

Cgj  - Coefficient  of  cospectrum  of  water  vapor  and  temperature 

e(T)  - Vapor  pressure  at  temperature  T (mb) 

. 2 

g - 9.8  m sec  force  of  gravity 
k - .35,  von  Karmann  constant 

L - Moni n-Obukho V length  (m) 
n - Optical  index  of  refraction 

P - Air  pressure  (mb) 

q - Vapor  density 

RH  - Relative  humidity  at  height  Z,  (%) 

Ri  - Bulk  Richardson  number 
R - .74 

T - Temperature  (K) 

- Temperature  at  height  Z 
T^  - Sea  surface  temperature 
Ty  - Virtual  temperature  (°K) 

Tva  “ Virtual  temperature  at  height  Z 
Ty^  - Virtual  temperature  at  sea  surface 
T*  - Friction  temperature 

U.  - Wind  speed  at  height  Z,  (usually  Z = 10  m)  (m  sec‘M 

d d 

U*  - Friction  speed 

Z.  - Height  (m) 

a 
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- Roughness  length  (m) 
a - Roughness  length  constant 
B - 4.7  constant  of  Barker  & Baxter  (1975) 
■ Heat  flux  profile  parameter 
Y^  - Momentum  flux  profile  parameter 
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APPENDIX  D 


2 

Humidity  Effects  in  within  the  Surface  Layer 

by 

H.  A.  Panofsky 


We  begin  by  defining  N* , by  analogy  to  T*, 


N 


U. 


(0.1  ) 


where  n'  represents  fluctuations  of  refractive  index  and  remaining 
terms  are  defined  in  the  body  of  the  report.  For  the  optical 
wavelength,  n'  is  related  to  temperature  T’  and  humidity  q' 
fluctuations  by  (Friehe,  1977) 


n ' = ( K.|  T ' + K^pq  ' ) x 


10 


-6 


(0.2) 


-1  3 

where  K,  = -1 .0°K  and  K,  = -.05667  cm  /microgram.  The  air 

. ^ - 3 - 3 

density*  p,  is  taken  to  be  1 .25  x 10  gm  cm 

Multiplying  (D.2)  by  w ' and  dividing  by  U*,  we  get 


N*  =-(!*+  71Q*)  X 10 
• - T.d  * 


-6 


(0.3) 


Next  we  introduce  Moni n-Obukho v scaling  with  the  scaling 

2 

function  N*  for  index  fluctuations.  Then,  proceeding  as  with  Cj 
we  get 


cl  = N*  G(z7L) 


(0.4) 


If  the  universal  functions  for  q and  T are  equal  to  each 
other  (as  they  seem  to  be  over  land)  G(z/L)  should  equal  g(z/L). 
This  can  be  seen  more  exactly  if  dissipation  of  refractive  index 
fluctuations  is  in  balance  with  the  production  of  these  fluctua- 
tions. Hence  we  finally  have 
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I 


C2  . (K3)2  c2  (1  . 

where  = 10  ® P/T^.  This  is  identical  to  the  result  of  Wesely 
(1976)  when  temperature  and  moisture  are  taken  to  be  perfectly 
correl ated . 


54 


